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This talk is based on Ref. [jl|]. The reader is referred to that paper for additional explanations. 

1. Introduction 

The energy between a static (i.e. infinitely heavy) quark and a static antiquark separated a 
distance r is known as the quantum chromodynamics (QCD) static energy, Eo(r), and is a basic 
object to understand the dynamics of the theory. As it is well known, one can distinguish a long- 
distance part and a short-distance part of the static energy. The long-distance part encodes the 
confining dynamics of the theory, whereas the short-distance part can be computed in perturbation 
theory. On the other hand, one can use lattice QCD simulations to compute the static energy in 
both (short- and long- distance) regimes. Here we want to focus only on the short-distance part 
of the static energy, where the perturbative weak-coupling approach is expected to be reliable. In 
particular, we want to compare state-of-the-art perturbative calculations with the lattice results. 
This comparison allows us to obtain a determination of the strong coupling a s , which is the main 
outcome of the present analysis. 

2. Perturbative calculation 

At short distances, Eq{t) is given at leading order by the Coulomb potential (with the adequate 
color factor), Eo(r) ~ —CpCCs/r (where Cf = (A£ — l)/(2N c ), and N c is the number of colors), 
and then we have corrections to this result. At present the static energy is known including terms 
up to order af +n ln n a s with n > [g g, |, |, g, 0, g] (a level of accuracy which we refer to as 
next-to-next-to-next-to leading logarithmic -N 3 LL-). The presence of lna $ terms in the expansion 
of the static energy is due to the virtual emissions of gluons with energy of order Eq (the so-called 
ultrasoft gluons), that can change the color state of the quark-antiquark pair from singlet to octet 
and vice versa []|, [Kj. Detailed expressions for the static energy at this level of accuracy were 
given in Ref. |g] (and references therein; see also [ |TT| ] for explicit numerical expressions), and we 
will not reproduce them here. 

3. Lattice computation 

The static energy has been recently calculated in 2 + 1 flavor lattice QCD [|l2|]. This com- 
putation used a combination of tree-level improved gauge action and highly-improved staggered 
quark action [|I3[]; it employed the physical value for the strange-quark mass m s and light-quark 
masses equal to m s /20 (corresponding to pion masses of about 160MeV). It was performed for 
a wide range of gauge couplings, 5.9 < j8 = 10/g 2 < 7.28. At each value of the gauge coupling 
one calculates the scale parameters ro and r\, defined in terms of the static energy Eq(t) as follows 



[14, 15] 



r^-U = 1.65, r^-U, = l. (3.1) 



The values of ro and r\ were given in Ref. [ 12] for each /3 . The above range of gauge couplings cor- 
responds to lattice spacings 1.909/ro < aT x < 6.991/ro. Using the most recent value ro = 0.468 ± 
0.004 fm ||T^] we get 0.805 GeV < aT x < 2.947 GeV. Thus we can study the static energy down to 
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distances r = 0. 14ro or r ~ 0.065 fm. For the comparison with perturbation theory the most relevant 
data set is the one that corresponds to j8 = 6.664,6.740,6.800,6.880,6.950,7.030,7.150,7.280, 
which is what we use here. The static energy can be calculated in units of ro or r\ . Since the static 
energy has an additive ultraviolet renormalization (self energy of the static sources) one needs to 
normalize the results calculated at different lattice spacings to a common value at a certain distance 
(or alternatively one can take a derivative and compute the force). The static energy is fixed, in 



units of ro, to 0.954 at r = ro (J12J. At distances comparable to the lattice spacing the static energy 
suffers from lattice artifacts. To correct for these artifacts we use tree level improvement. That is, 
from the lattice Coulomb potential 

■0zD oo (k o = O,k)e itr , (3.2) 

we can define the improved distance rj = {47iCi{r))~ x for each separation r. Here Dqq is the tree 
level gluon propagator for the a 2 improved gauge action. The tree level improvement amounts to 



replacing r by rj []16J. Alternatively following Ref. []15|,|17|] we fit the lattice data at short distances 
to the form const — a/r + or + a'{\/r—\/rj) and subtract the last term from the lattice data. Since 
the data at the shortest distances that we use (for each j3) correspond to a separation of one lattice 
spacing, it is important to check that the way we are using to correct lattice artifacts is working 
properly. In that sense, we have found that both methods of correcting for lattice artifacts lead to 
the same results within errors of the calculations. Furthermore, the static energies calculated for 
different lattice spacings agree well with each other after the removal of lattice artifacts, and when 
one puts all the data together it seems to lie on a single curve, even at short distances, indicating 
that the above procedure of removing the lattice artifacts works. 

4. Comparing lattice and perturbation theory: a s extraction 

We can now compare the lattice results with the perturbative expressions, and use the com- 
parison to extract the value of the QCD scale A^r (in the MS scheme). In order to obtain this 
extraction, we assume that perturbation theory (after implementing a cancellation of the leading 
renormalon singularity) is enough to describe lattice data in the range of distances we are consid- 
ering (we use lattice data for r < 0.5ro, and since we have lattice data points down to r = 0.14ro, 
this means that we are studying the static energy in the 0.065 fm< r <0.234 fm distance range, in 
physical units). Then we search for the values of Aj^t that are allowed by lattice data; the guiding 
principle to do that is that the agreement with lattice should improve when the perturbative order 
of the calculation is increased. 

As it was already mentioned above, in the perturbative calculation of the energy one needs to 



implement a scheme that cancels the leading renormalon singularity [18, 19]. This kind of schemes 



introduce an additional dimensional scale in the problem (that we denote as p). We implement the 



renormalon cancellation according to the RS -scheme described in Ref. [20]. The static energy in 
this scheme is given by 

£ RS (r,p) = -RSsubtr.(p), (4.1) 

where the subtraction term on the right-hand side cancels the leading renormalon singularity of 
E^ s (r); the explicit expression for RSsubtr.(p) is given, for instance, in Eq. (7) of Ref. [ lip. The 
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scale p has, in this case, a natural value which corresponds to the center of the range where we com- 
pare with lattice data (i.e. around 1.5 GeV), but any value around that one cancels the renormalon 
and is, therefore, allowed. 

4.1 Central value for foA^ 

To obtain our central value for foA^ we use the following procedure: 

1. We let p vary by ±25% around its natural value. 

2. For each value of p and at each order in the perturbative expansion of the static energy, we 
perform a fit to the lattice data (roA^ is the parameter of each of the fits). 

3. We select those p values for which the reduced % 2 of the fits decreases when increasing the 
number of loops of the perturbative calculation. 

Then we consider the set of roA^g values in the p range we have obtained and take their average, 
using the inverse reduced x 2 of each fit as weight. From that, we obtain our central value for 
r o^MS- The value we obtain at 3 loop with leading resummation of the ultrasoft logarithms is 
r o y ^Ms = 0-70, which will be our final number for the central value. The perturbative expressions 
at N 3 LL accuracy (i.e. 3 loop with sub-leading resummation of the ultrasoft logarithms) are also 
known (as mentioned before) but in this case an additional constant appears in the expressions (due 
to the structure of the renormalization group equations at this order This additional constant 
would also need to be fitted to the lattice data (i.e. one has a two-parameter fit in this case). When 
we do that we find that the % 2 as a function foA^r is very flat, and we cannot improve the extraction 
by including these higher order terms. In principle, more precise lattice data, and/or data at shorter 
distances might allow for an improvement in that respect. 

4.2 Error estimate 

Having obtained our central value for foA^r we now need to assign an error to it. We want the 
error to reflect the uncertainties associated to the neglected higher-order terms in the perturbative 
expansion of Eo(r). To achieve that, we consider: (i) the weighted standard deviation in the set 
of p values we found above, and (ii) the difference with the weighted average computed at the 
previous perturbative order. (Note that, starting at two-loop order, one can decide whether one 
wants to perform the resummation of the ultrasoft logarithms or not. To assign the error we take 
whichever difference -with or without resummation in the previous order- is larger. This amounts 
to not making any assumption about the necessity or not to resum these logarithms). We then add 
the two errors linearly (term (ii) turns out to be the dominant one). 

Additionally, we also redo the analysis with alternative weight assignments (/rvalue, and con- 
stant weights); we obtain compatible results. In the final result, we quote and error that covers 
the whole range spanned by the three analyses. As an additional cross-check, we can compare the 
analysis performed with the static energy normalized in units of ro (our default choice) and the 
one with the static energy normalized in units of r\ . We find that the two analyses give consistent 
results. (Note that the values for the static energy in both cases -i.e. in units of ro or r\ - come from 
the same lattice data set in terms of r/a; but the error analysis in the normalization of the energy 
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for each lattice spacing is different in the two cases. Therefore, one cannot obtain Eo(r/r\) from 
Eo(r/ro) by a trivial rescaling, and it is in this sense that the analysis with the scale r\ provides a 
cross-check of the result). 

4.3 Final result for a s 

Our final result reads 

r AMs= 0.70 ±0.07, (4.2) 

which corresponds to 

a s {p = 1.5GeV,n/ = 3) =0.326 ±0.019 a s {M z ,n f = 5) = 0.1156+g;^, (4.3) 

where we used the ro value from Ref. [12] to obtain a s (1.5GeV) from Eq. (|^2|), and then evolved 



it to the Z-mass scale, Mz, using the Mathematica package RunDec [ |21[ ] (4 loop running, with 



the charm quark mass equal to 1.6 GeV and the bottom quark mass equal to 4.7 GeV). 



4.4 Comparison with other recent a s determinations 

There are several recent determinations of a s that also employ comparisons with lattice data. 
These include analyses that use: observables related to Wilson loops (but not the static energy) 



[22, E3p, moments of heavy quark correlators [24, 23], the vacuum polarization function [ESQ, the 



Schrodinger functional scheme Q26Q, and the ghost-gluon coupling Q27|[. They deliver numbers that 
are mostly compatible with our result, although our central value is a bit lower than those of the 
other lattice determinations (see Fig. |l| for a graphical comparison). 

In Fig. || we compare our result for a s (Mz) in Eq. (43) with a few other recent a s determi- 
nations that use other techniques (i.e. non-lattice determinations); we include results coming from 
T decays, thrust, and partem distribution function (PDF) fits 1 , along with the Particle Data Group 
(PDG) average (the comparison is not exhaustive, we just show a few other recent results in the 
figure, meant to illustrate where our result lay with respect to recent a s extractions). 

It is also worth remarking that our a s determination is performed at a scale of around 1.5 GeV, 
and therefore constitutes an important new ingredient to further test the running of the strong cou- 
pling (see Fig. 2 of Ref. [[[]] for a graphical comparison of different determinations of a s as a 
function of the energy scale Q where they are performed). 



5. Conclusions 



To summarize, in this work we have compared perturbative calculations for the QCD static 
energy at short distances with lattice computations. We find that perturbation theory (after can- 
celing the leading renormalon singularity) is able to describe the short-distance part of the static 
energy computed in 2 + 1 flavor lattice QCD (see Fig. 1 of Ref. ^ for a comparison of the different 
orders of accuracy of the perturbative result and the lattice data). We exploited this fact to obtain 

1 Note that the errors in the results from PDF fits do not include effects from the unknown higher-order perturbative 
corrections. This theoretical uncertainty is difficult to assess and has not been addressed in detail so far. It is expected to 
be roughly of the same order as the quoted errors. 
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Figure 1: Comparison of our result for a s (Mz) with other recent lattice determinations. The references are: 
HPQCD @, JLQCD [||], PACS-CS [§j|, ETM @- 
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Figure 2: Comparison of our result for a s (Mz) with a few other recent a s determinations. We include 
results from z decays (Boito et al. [|||; Abbas et al. [Q; Pich @), thrust (Abbate et al. [||]), and PDF 
fits (ABM 11 ||], MSTW [fH, NNPDF @), along with the PDG average ^ for reference. 
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a determination of the strong coupling a s . Our extraction is at three-loop accuracy (including re- 
summation of the leading ultrasoft logarithms) and is performed at a scale of 1 .5 GeV. When we 
evolve the result to the scale M z it corresponds to a s (M z ) = 0. 1156^0 0022- 
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